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SINDA/FLUINT has been found 


abstract 

nd to be a versatile code for modeling aerospace systems involying 
sing ,e flowed 5f -ode. of ^ “ht 

Johnson Space Center (JSC) for modeling , con( j e nsation two-phase flow pressure drop, 

have been utilized xnclude tr^sient simulTO ^ been u$ed at JSC to supp ort the 

^Keral o, these applications of 

SINDA/FLUINT at JSC are described in this paper. 

SINDA/FLUINT is being used extensively t0 

phase ammonia loops of the Space Station e U bmodelsto very detailed subsystem models. 

between the multiple fluid loops which comprise the Space Station ATCS. 

JSC has also developed several detailed models of ^ SC* 

Space Station ATCS. One model, which has approximate y^ tw J fl u jd loops connected only 

provides a representation of the ATCS low tempera models are used to determine parameters such 

SjSSSf ASSESS and total heat rejechon for the baseline 

, design as well as proposed alternate designs. 

SINDA/FLUINT has also been used at JSC as a design 
gasses. One model examined the pressurization and internal cooling. Another 

variety of operating conditions including convection man euverine units This model included 

model predicted the performance of a new generate _of SJJLr equations. The results of 

high pressure gas depressunzation, internal eas bottles and also to point to areas 

both models were used to size components, such as the heaters and gas ootucs «u.u 

where hardware testing was needed. 
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INTRODUCTION 


Use of the Systems Integrated Numerical Differencing Analyzer (SINDA) analysis tool has 
expanded steadily over the years since its origin in the 1960's. The Fluid Integrator (FLUINT) code 
added significant fluid system analysis capabilities under a NASA Johnson Space Center (JSC) contract in 
1985 (ref. 1). Since then, analysis applications for the code have increased in scope as well as in number. 
At the Johnson Space Center, SINDA/FLUINT has been used to solve both steady state and transient 
thermal/ hydraulic problems involving single and two-phase fluid flow with heat transfer. Conduction, 
convection and radiation have been modeled in simple and complex systems alike. Models have ranged 
from a few nodes to a few thousand nodes. SINDA/FLUINT has been used at JSC to support the Space 
Shutde, Space Station Freedom (SSF) and other advanced technology programs. The code has been used 
in the conceptual design, detailed design, test and performance verification phases of these programs. 

The Space Station Thermal Control System (TCS) consists of both passive thermal control features 
and active fluid loops to transport heat The power system uses a single-phase ammonia cooling system 
and the habitable modules use single-phase water loops to gather heat from racks filled with equipment 
The internal thermal control system (TTCS) water loops transfer their heat to one of three external thermal 
control system (ETCS) two-phase ammonia loops which dump the heat to space. Figure 1 shows a 
schematic representation of the Permanently Manned Capability (PMC) Space Station Freedom and 
indicates the location of some of the TCS features. Thousands of SINDA/FLUINT CPU hours have been 
logged analyzing the various Space Station Freedom Thermal Control Systems during the design phase, 
and this use will continue into the operational phase since SINDA/FLUINT is the official code used for 
Space Station Integrated Thermal Analysis. SINDA/FLUINT submodels developed by various system 
designers are integrated and run at JSC in order to analyze integrated system performance under various 
nominal and off-nominal conditions. 

Other more detailed SINDA/FLUINT models have been built of particular components of the SSF 
Active Thermal Control System (ATCS), particularly the ETCS radiators. These radiators have undergone 
several redesigns in order to improve performance and reduce weight, cost and assembly time. 
SINDA/FLUINT has been used in these design trade studies in order to evaluate potential alternate 
designs. One such study, which will be described here, looked at ways to prevent the radiators from 
freezing when subjected to an environment temperature substantially below the freezing point of the 
ammonia working fluid. 





Figure 1: PMC Configuration of Space Station Freedom 
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SPACE STATION FREEDOM INTEGRATED THERMAL ANALYSIS 

SINDA/FLUINT has been utilized in 
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conditions include nominal and off-notninal operations. 

Analysis of nominal operations generally analysis 

^^Sl'oS^ve ‘S^A/FLUINT models will be briefly described and examples 

of problems solved will be given. 

SINDA/FLUINT Model Descriptions 

The SINDA/FLUINT models utilized 

habitable modules (internal thermal control system ' - ^dth rmod 1 1 ™J nodds rf ^ intcrnal 

ammonia heat transport system (external thermal control syste habitation and laboratory modules, 

both adetailed Sid a simplified version of the models have been developed. 

The modeling philosophy utilized for die S(SS 
the exit temperature of each experiment rack °ff^ P , modeled. The transient response time is also 
conceivable configurations and transients .can be y ^ etCS. The models do not consider the 
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nodes. 

‘^hrXrul ^u»cSTti^ng cons, tiered by Che 

International Partners. 
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The U.S. Hah/Lab system uses a System Flow Control Assembly tq maintain a constant pressure 
drop between the rack supply and the rack return lines. The maintenance of the constant "delta-p" helps 
ensure adequate flow to each rack. The flow to a particular rack is controlled by the Rack Flow Control 
Assemblies which monitor either temperature or flowrate, depending on the particular type of control 
desired. Variation of a rack's flowrate is necessary to conserve power since not all racks will require 
cooling at all times. The U.S. Hab/Lab system also includes a regenerative heat exchanger that is used to 
ensure that the water entering the MT branch is above 60°F (15.5®C). The 60°F setpoint was chosen since 
the Space Station maximum allowable dewpoint is 60°F, thus precluding condensation on the MT branch 
of plumbing and saving the expense of insulating the entire system. Each Interface Heat Exchanger (IHX) 
has a bypass line that is used to maintain a desired overall IHX outlet temperature. 



Figure 2. United States Habitation/Laboratory Functional Schematic 


The FLUINT models of the habitation and laboratory modules are similar, the primary difference 
being the number of payload and system racks. The level of detail is such that all racks are modeled as 
well as important control systems. Figure 3 provides a schematic showing the number and location of the 
FLUINT components for the U.S. Laboratory module. 

For the U.S. Lab model, most fluid line lengths and sizes have been defined and have been 
incorporated. The locations of elbows and some connections have not yet been defined. The control valve 
characteristics, such as k-loss and flow areas, have been defined for the Rack/System Flow Control 
valves. The characteristics for the remainder of the valves are not yet defined. The Rack/System Flow 
Control valve is a sculpted ball valve with a variable k-loss. This action is modeled within FLUINT by 
altering both the valve angle and the k-loss value. The thermal mass of the equipment in each rack is 
accounted for by including a SINDA thermal node equivalent to 30 lbm (13.62 kg) of stainless steel, 
though the actual mass of equipment in each rack is currently unknown. 
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The model of the U.S. Hab maintains the same level of detail, though the number of FLUINT 
components is reduced since the U.S. Hab has fewer racks. 

Resource Node System and Model D escription „ 

The resource node contains the primary command and control avionics tor the Space Station, i ne 
resource node ITCS cools the resource node 2, resource node 1, airlock and Pressurized Logistics Module 
(PLM). The PLM can be located at either resource node 1 or resource node 2. Since the airlock, resource 
node 1 and PLM are attached at different stages of the Space Station construction, and possibly at different 
locations for the PLM, several distinct configurations of the resource node ITCS exist This variability 
forces special modeling considerations for the resource node and will discussed below. The architecture 
of the resource node is the most developed at this time since the resource node will be the first habitable 
module launched. Therefore, the FLUINT model of the resource node is the most complex and also the 
most useful for predicting ITCS performance. 

The ITCS for the resource node is a simplified variant of the ITCS for the U.S. laboratory module 
discussed above. Unlike the U.S. Hab/Lab systems, the equipment in the resource node is normally 
always on. Therefore, the need to provide a variable flowrate is not present. The ITCS developers have 
chosen to use an orifice to control the flow to each rack with the coolant pumps set to a constant speed. 
There are some active components associated with the resource node, the primary ones being rack flow 
control valves associated with the PLM. These valves will modulate to maintain a constant exit 
temperature for each PLM rack. Also associated with the PLM is a modulating delta-pressure valve which 
acts to maintain a constant pressure for the inlet header to the PLM racks. 

The detailed FLUINT model of the resource node uses the same component descriptions and 
modeling as the U.S. Hab/Lab due to common hardware. Such items as IHX bypass valves, RHX 
modeling, and pump modeling are identical. The modeling of the rack is also the same. Modeling specific 
to the resource node involves the level of detail of the plumbing and the variability in the configuration. 

The geometry of the resource node is established sufficiently such that three-dimensional isometric 
drawings, with dimensions, have been generated. These drawings allow for a more accurate definition of 
ITCS line length. Also, the pressure drop effect of elbows and bends can be estimated using the effective 
length approach. Though the FLUINT model does not yet consider the effect of fluid connectors, such as 
tees and quick disconnects, the actual connectivity of the ITCS can now be included into the model, with 
the pressure drop data included when available. 

The variability of the resource node configuration results in four unique configurations, not 
counting the position of the PLM. The multiple configurations of the ITCS are accounted for by utilizing 
the capability of FLUINT to include sections of input from predefined external files. Each module is 
contained in a separate input file and included as necessary to create the ITCS connectivity of the particular 
configuration. The use of external files allows for easy maintenance of the various portions of the 
resource node ITCS model. Should a particular module change, only the model of that module need be 
altered to make the change effective for all configurations. The alternative would be to make changes to as 
many as eight separate models, which could lead to errors in implementation. 

International Partners System and Model Description . 

The International Partners, the European Space Agency (ESA) and the National Space 
Development Agency of Japan (NASDA), are currently in die process of re-designing their respective 
ITCSs. The baseline design for both is a one-loop system with temperature controlling bypasses. The re- 
design efforts are due to problems in meeting Space Statical thermal load management goals. The final 
configuration will not be known for several months, however some of the basic components will probably 
still exist. 

The ESA ITCS is similar in concept to the resource node ITCS. The coolant flow to each primary 
payload rack is controlled via an orifice. In parallel with the racks is an avionics heat exchanger. The 
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inputs. 
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phase coldplates, and associated plumbing ■ ® 0 _ , OTC Moderate Temperature (MT) loop 

SS M.?Q r^iffifu^onal scheme for ,he Man 

Tended Configuration of SSF. 



Figure 4. ETCS Functional Schematic for the Man Tended Configuration 
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The system operation consists of the RFMD suppling slightly subcooled ammonia liquid to a 
distributed system of heat exchangers. The flow to each heat exchanger, and therefore the ultimate heat 
acquisition capability of each heat exchanger, is controlled by a cavitating venturi placed at the inlet of each 
heat exchanger. The downstream pressure is controlled so that the ammonia will be very close to the 
saturation point after exiting the cavitating venturi. With the addition of heat from the internal modules, the 
ammonia begins to boil, enhancing the heat transfer process. The heat transferred to the ammonia is 
controlled so that the maximum ammonia vapor quality is no greater than 90%. The resulting two-phase 
mixture returns to the RFMD. 

The RFMD consists of a stationary outer housing with a rotating inner housing that acts to separate 
the two-phase mixture into vapor and liquid components (Figure 5). The ammonia vapor is plumbed to 
the condensing radiators where the acquired heat is rejected to space. The condensate then returns to the 
RFMD where its temperature is increased and then pumped back out to the module's heat exchangers. The 
pumping action of the RFMD is via immersed pitot probes which convert the rotational energy into a static 
pressure rise, thus providing the necessary pumping power. The RFMD maintains an internal liquid level 
by means of a level pitot which acts to pump excess fluid to the bellows accumulator. 
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Figure 5. Rotary Fluid Management Device Schematic 


Since the ETCS is a two-phase system, the operating pressure controls the operating temperature. 
The Back Pressure Regulating Valve is designed to passively control the pressure in the RFMD. The 
BPRV does this through a combination of spring and servo action (Figure 6). The system is also capable 
of changing operating temperatures by altering die relative spring forces inside the BPRV, thus temporarily 
affecting the heat rejection rate and raising or lowering the system pressure. 


The modeling of the ETCS can be broken down into two distinct portions; the first being the 
modeling of the pump module assembly (RFMD, BPRV, and bellows accumulator) and radiator and the 
second being the remainder of the system (main plumbing, cavitating venturies, and module heat 
exchangers). The model of the pump module assembly consists of a detailed description of the RFMD, 
BPRV, and bellows accumulator. The model considers the physics of the RFMD by determining 
parameters such as the induced gravity head due to the rotation, pumping power from the pitot probes, and 
internal liquid levels. The model of the BPRV uses the internal geometry, including spring forces and 
internal bellows areas, to calculate the forces on the primary pressure control valve. TTte condensing 
radiator model is of moderate detail and accounts for thermal interactions between the parallel flow 
passages. 
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Figure 6. Back Pressure Regulating Valve Schematic 
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Summary of Analyses Performed with the SSF Integrated ATCS Models 
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the ETCS which could have caused an overload of the ETCS houowng a re 6^ Node Also ^ 

SJM problem on a portion of the plumbing 
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that was un-insulated. The models were used to evaluate various options and to quantify the heat load 
value at which condensation became a potential problem. 



□ TWO-PHASE AMMONIA 
0 HEAT ACQUISTION POINT 

Figure 7. FLUINT Nodalization of the Moderate Temperature Loop for PMC 


For the International Partners, the models identified a similar mal-distribution of module waste heat 
that violated Space Station operating specification. Several IP proposed options were modeled and 
evaluated as to the option's ability to meet station requirements. Additionally, the models have been used 
to evaluate control system stability. 

In the integrated mode, the models have identified several key system issues and been used to 
analyze system operation at off-nominal conditions. One system issue identified was the inadequate heat 
load placed on one of the loops of the ETCS during early operations. The ETCS nominally consists of 
three loops - two Low Temperature and one Moderate Temperature. The Low Temperature loops are the 
first activated during the early stages of space station construction, when the total heat loads are the lowest. 
The pump module assembly requires a minimum amount of heat load (about 1.5 kW) to ensure adequate 
operation of the BPRV. Integrated analysis showed that only about 1 kW of thermal load was present on 
one of the Low Temperature loops during certain portions of station construction. The integrated models 
woe used to devise methods to shift heat load from one Low Temperature loop to the other by altering the 
operating temperatures. Also identified by analysis in the integrated mode was the possibility of ammonia 
freezing in the condensing radiators due to low heat loads and cold external environments. The radiator 
manufacturer has confirmed this possibility and the design is being modified to account for freezing. 

Off-nominal analysis has involved investigating the impacts to the system when operated at 
elevated temperatures. Operation at elevated temperatures enhances the heat rejection capability but can 
adversely impact the atmospheric temperatures inside the habitable modules. Since this change impacted 
the habitable modules, a coordinated analysis effort with the module providers was necessary. The ease 
with which the SINDA/FLUINT models can be changed allowed for rapid evaluation of proposed 
operating points. 
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Figure 8 : Detail of the Tube Extrusion 
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Two separate sets of manifold tubes (manifolds A and B) carry vapor to and retum condensate from each 
panel. Alternate tubes in each panel are connected to one set of manifold tubes. A diagram of the flow 
distribution to each panel in each ORU is shown in Figure 9. The manifold tubes at each end of the 
radiator panels are enclosed by manifold cover plates made of 6061 -T6 aluminum and are 5.28" (13.4 cm) 
wide. The overall length of the panels including the manifold covets is 131.46" (333.9 cm). In order to 
provide support for the face sheets and increase the strength of the panels, the space between face sheets is 
filled with a layer of hexcel honeycomb material with a density of 3.1 lb/ft^ (49.6 kg/nw). The 
honeycomb is made of 5052 aluminum foil 0.0007" (0.018 mm) thick. 

A moderate temperature (MT) ammonia loop (62° F, 16.6° C) flows through all port side radiator 
tubes. Two low temperature (LT) ammonia loops (35° F, 1.66° C) flow through alternating tubes in the 
starboard side radiators. No mass transfer occurs between die two LT ammonia loops, but they may 
communicate thermally by conduction through the face sheet 

The size of the plumbing of SSF has been optimized to minimize the weight while maintaining die 
pressure drop in lines at an acceptable level. Therefore, different size tubes are used in the LT and MT 
loops. Table 1 shows a summary of the plumbing sizes used in this study. The 3 ORU's at each side are 
pre-integrated on a section of the SSF truss in a folded position for easy transportation and can be 
deployed on orbit. This requires that the radiator panels in each ORU be connected by flexible tubing that 
can tolerate folding and unfolding. The manifold tubes of each panel in each ORU are connected to the 



ORU 1 ORU 2 ORU 3 

Figure 9 : Schematic of Manifolds and Radiator Panel Layout for one Wing 
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Table 1: Plumbing Size Summary ' 
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coldplates. The estimated pressure drop associated with each WXor^r was mooc cu y 
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Figure 10 : Arrangement of the Manifold Tubes and Flex Hoses 


The rotary fluid management device (RFMD) was modeled by two plena - one representing the 
warm end set at saturation pressure corresponding to the desired set point temperature, and the other 
representing the cold end set at subcooled condition. A pressure difference equal to the RFMD end-to-end 
pressure difference (as reported by the manufacturer) was imposed on two plena. A saturated liquid flow 
rate equal to the evaporator flow rate plus the RFMD bearing flow and the back pressure regulating valve 
(BPRV) servo flow was extracted from the warm end plenum. The vapor flow rate out of the RFMD was 
calculated dynamically by performing an energy and flow balance on the RFMD. The vapor flow rate is a 
function of the instantaneous heat load, set point of the RFMD, evaporator and condenser return flow 
conditions, and bearing and servo flow rates. The vapor flow mixes with the BPRV servo flow before 
entering the radiator panels. The servo flow rate was assumed to vary linearly with total heat load on the 
system from 0.022 GPM at maximum heat load to 0.015 GPM at minimum heat load (3.35 kW). The 
BPRV is designed to maintain the RFMD warm end pressure at the set point saturation pressure by 
regulating the vapor flow rate. The BPRV opening (which determines BPRV pressure drop) is then 
adjusted to allow sufficient vapor flow out of the RFMD. The model of the BPRV performs the same 
functions for the normal operations of the system. The BPRV pressure drop in the model varies with 
vapor flow rate in order to balance the RFMD end-to-end pressure drop. The model will signal if the 
BPRV reaches its maximum opening and can no longer control die set point. 

Several versions of the radiator model with different levels of detail were developed for different 
analyses. The alternating flow tubes in each radiator panel connect either to LT loop A or B. The flow 
tubes and other lines and components corresponding to each loop were represented in a separate submodel 
since no flow mixing occurs between the two loops. Each panel was then represented by models of one 
flow tube from loop A and one from loop B which were thermally connected through conduction in the 
face sheet nodes. Each flow tube was modeled by a HX macro with 10 segments in the direction of the 
flow. The ammonia in each section was represented by a JUNC and the flow tube section was represented 
by a TUBE with each lump downstream of each tube section. 
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Condensation Heat Tra nsfer Correlat ion , Rohsenow's conclation (ref. 2). 
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Thermal Network ‘‘ . 

The radiator panel face sheets and extrusion tubes were divided into 10 strips along the panel 
length. The face sheet area corresponding to each tube segment was 4.73" (12.0 cm) wide, 12.09" (30.7 
cm) long, and centered above the tube. The distance of 4.73" (12.0 cm) on the face sheet between the two 
adjacent tubes modeled was divided into 9 nodes. This nodal break down allowed for determination of the 
temperature profile of the face sheet between two tubes and eliminated the need for a pre-determined fin 
efficiency for the radiator face sheet. The nodal breakdown of the radiator face sheet is shown in Figure 
12. The mass of each tube segment was represented in the capacitance of die tube wall nodes which was 
thermally tied to the fluid lump in that tube segment via a heat transfer tie. The flow tube inner wall node 
was connected to the extrusion inner surface via conductance through the stainless-steel and tube adhesive. 
The overall conductance of the extrusion tube, from its inner surface to its interface surface, was calculated 
from a separate detail SIND A model of the extrusion tube. The extrusion interface node was connected to 
die face sheet node via conduction through the adhesive. The extrusion tube nodes were also connected 
axially to allow for axial conduction in the subcooled part of the tube. The face sheet nodes were 
connected both laterally and axially. The thermal network of one segment of two adjacent tubes is shown 
in Figure 13. 

Manifold cover nodes 



tube spacing 4.73" 


Figure 12: Nodal Breakdown of Radiator Facesheet 


236 
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flow in alternating tubes of the radiators in the same considered in which the flow 

tolerance of the radiators, an iM c ^«Sof Ae^dfolds for the parallel . 

L ■* - te presCTted ta 

addition to the results for the baseline design. 
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condensate 
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manifold A 
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manifold B 


Vapor 
manifold A 


Vapor 
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Parallel-flow 


vapor 


Liquid Vapor 

manifold A manifold B 


Vapor Liquid 

manifold A manifold B 



Counter-flow 

Figure 14: Schematic of the Radiator Panel Manifolds for Parallel and Counterflow 


Detailed Radiator Model Analysis Results 

The models were analyzed at steady-state and transient conditions for a variety of cases as 
described below. The Space Station configuration analyzed was the Man Tended Capability (MTC) stage 
with only one of the radiator ORUs deployed. 
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^^^fficconhlS &de. was analysed at 

temperature of -50°F (-45.5°C). The heat lo ^ d J^?j C 5 u xheresults indicated that 2-phase ammonia with a 
the mo loops wwUentieg $ « *g™ liquid a, -32£F 

vapor quality of 0.967 at 29.3 F (-1.5 C) condensing and the rest was subcooling the flow. 

(-35 9°C). 72.6% of the length of each flow tube was conaensg overall fin efficiency is the 

equal to the temperature of the hottest face sheet node. 


Table 2 : Results of the MTC LT loop analysis 
at -50° F sink temperature 



Flow Rate 
Lb/hr 

Hut Load 
CKW] 

Inlet 
Temp [F] 

Inlet 

quality 

Outlet 

r«*pm 

Condensing 

Fraction 

Loop A 

0.498 

7.5 

29.37 

0.967 

-32.7 

0.726 

Loop B 

0.498 

7.5 

29.37 

0.967 

-32.9 

0.726 


=7x^3^ - "" “ - 839 ” F 

(-64.4°C). 

The model was run for the same cold case conditas ; ^“SSotlde? * 
configuration to examine the effectiveness of jhc ^perature was 84.5"F <47=0 
temperature. The results shown in Table 3 indicated flow tube for the parallel and 

Sfor the counter-flow case. The ammom ‘ a ™ u 
counter-flow cases are shown in Figure 15 . The re^ ts ind r __ densate f^ w at the outlet was much 
completed faster in the counter-flow case, adjacent tube The coldest ammonia temperature in 
less Sue to heat leak ^ the 2_Phase Sri&etfthe panel. The panel overall fin 

the counter-flow case was -60 F (-5 1 . l y j rt cc rpcnectivelv 

effidencies for the parallel and counter-flow cases were 0.54 and 0.55, respectively. 


Table 3 : Results of the MTC LT loop analysis 
at -95° F sink temperature 



Flow Rate 
Lb/hr 

Hut Load 
[KW] 

Inlet 

TemptF] 

Inlet 

quality 

Outlet 

TenipIF: 

Condensing 

Fraction 

Parallel 

flow 

0.457 

7.5 

29.1 

0.97 

-83.9 

0.48 

Counter- 

flow 

0.528 

7.5 

29.1 

0.97 ' 

0.6 

0.29 
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Figure 15: Ammonia Temperature Profile in Flow Tube A for Parallel and Counterflow 


Radiator Freeze Prevention Studies 

Under low heat loads and very cold environments when the effective sink temperature falls below 
the freezing temperature of ammonia (-108°F, -78°C), the fluid in the radiators may freeze. This leads to 
very high local pressures inside the flow tubes when frozen ammonia begin to thaw. During the SSF 
ETCS design process, the problem of radiator freezing (and subsequent thawing) was addressed. One 
approach was to try to increase the radiator condensate return temperature in cold cases without limiting its 
heat rejection capability in hot cases. Other approaches allowed the ammonia to freeze and thaw under 
controlled conditions. Several of the options which were investigated using the detailed SINDA/FLUINT 
radiator models are presented here. 

Freeze Prevention bv Radiator Pointing 

The MTC configuration of SSF was analyzed with two heat loads at transient conditions using a 
cold case environment (orbit angle Beta=0°, beginning of life surface properties) and different radiator 
orientations. The results are shown in Table 4. The coldest environment was the "edge-to-sun on the sun 
sideY'edge-to-earth on the dark side" (ETS/ETE) orientation, which resulted in radiator freezing for some 


Table 4 : Summary of the results of MTC configuration transient analysis under 
cold environment at [3=0, for different radiator panels orientation 


r PANEL 

ORIENTATION 
l sunside/darkside , 

NO. of 
Panels 

'jH it m il 1 1 1 il i ^ 

Condensate return temp. 

['FI 

. 44) kw 1 11.2 kw A 

ETS/ETE 

8 

freezing 

freezing 

ETS/FTE 

8 

-106 

-95 to -100 

FTE/FTE 

8 

-29 to -51 

- 

45°-to-earth 

8 

-62 to -75 

-33 to -81 
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parts of the orbit for both heat loads. The "face-to-earth on 

SSsSssasaMssas^^aa- 

woe much colder. Therefore, this option alone was not sufficient to prevent raoiator 8 

Frn~^ PrcycTl don bytn a^si^^v^^cwducte d ^S^ nine the feasibiUtytf other 

options such^as increasing the BPRV liquid' servo feu^lSd 

to prevent rediator freezing (ref. 4). n,e.4^t oold ^ /'X^d,tK^Sto?tacreasing 
TEA=0°. shown in Figure 16, was used in the study. The resulB^estu^^ trai^^ ^ 

BPRV servo flow had nosignmcrneffecKM^OTMdM Staeasedb^ tttoes its normal value. 

^ S^dtoaied 

(-67.8°C). Use of that much power would have an unacceptable impact on S . 



Freeze Prevention with Counts r-flpw J&faga i s . another option that was considered since it 

The counter-flow radiator design (described above) was anotner opuon u«oi 

does not increase the weight of the radiators and does not require any additi°n^ c^I»^nts ^ °P 

S“t way to alsF design ""^“^w^mluto 
freeze and thaw without bursting. For the SSF case, the rreeze/tn gn annlications the 

works even when one of the low temperature loops is not PP 

counter-flow option may be the best design solution to deal with radiator freezing. 
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A case was run using the MTC model at low loads (3.35 kW per loop), using the coldest panel 
environment for Beta=0 and TEA=-45° as shown in Figure 17. The results of the analysis using the 
parallel-flow (baseline) option are shown in Figure 18. The results showed that freezing of ammonia 
occurred at sub-freezing sink temperature parts of the orbit The results of the analysis with the counter- 
flow model are shown in Figure 19. Unlike the parallel flow case, the coldest fluid temperature in the 
counter-flow design occurs at the middle of the radiator panels since the condensate return flow at the 
outlet of each tube is warmed up by the 2-phase flow at the inlet of the adjacent tubes. Figure 19 shows 
the radiator condensate return temperature as well as the radiator coldest fluid temperature and indicates 
that the freezing problem was eliminated by the counter-flow option. 


|3 = 0°, o/e=0.15/0.92, TEA = -45° 




Time [Hour] 


Figure 18: Parallel-flow Radiator Exit Temp. (Beta=0°, TEA=-45°, 3.35 kW per loop) 
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Figure 19: Counternow Radiator Temperatures (Beta=0°, TEA=-45°, 3.35 kW per loop) 


Detailed Radiator Models Summary 

SIND A/FLU INT is a powerful tcx 
The detailed models developed at JSC for 
essential in predicting performance of the 
practical interest. The models have been 
environment transients and radiator freezing. 


for detailed analysis of complicated thermal and fluid systems, 
he analysis of the Space Station ETCS radiators have been 
adiators for many nominal and off-nominal conditions of 
__ j /Ucton nmhlems as response to oroitai 


IDEAL GAS SYSTEMS 

From thermodynamics, when a gas is compressed or expanded isentropicahyand atUabamany, its 
temperature will increase or decrease, as described by the integrated, isemropi ■ , j ^ j aw 

^S^comprcssion) or the first law of *ermod = cs ^Is^il? 

SSSS sSul^ TJSSSSm must be developed. The SINDA/FLUINT progtam out be 
used to solve these problems. 


The Basic SINDA/FLUINT Model 


Figure 20 presents the schematic of a basic SINDA/FLUINT mo*l Ite . can 1 £ 


removed 

ed 

i convection 
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heat transfer (HTU) tie. The value of this tie must be determined for each specific application. System 
flow rates are controlled by the MFRSET (Mass Flow Rate SET). The plenum (PLEN) is used to 
represent an infinite source of supply gas or a reservoir to dump the exhaust gas. The TANK is connected 
to the PLEN by an STUBE (Short TUBE) which can, if User Logic is employed, represent a pressure 
drop device (such as a regulator or an orifice). The working fluid is modeled as an ideal gas by using 
8000 series fluid data blocks. 

To ensure that this model was properly developed, the predicted final temperatures were compared 
against the theoretical results for adiabatic filling and isentropic, emptying processes for dry air. For a 
tank emptying process with constant thermophysical properties, the equation for die final temperature is 
(ref. 5), 



where Tf is the final temperature, Tj is the initial temperature, Pf is the final pressure. Pi is the initial 
pressure and k is the ratio of the specific heats. 

For a filling process, the following equation may be used (ref. 6), 


kPfTinTj 

lf "(Pf-Pi)Ti + kPiT in 


[ 2 ] 


where Tin is the temperature of the inlet gas. 

For the range of parameters to be considered here, the predicted results from the SINDA/FLUINT 
model were within 0.5 °F (0.3 K)of the theoretical results of both equations. The slight discrepancy is due 
to use of constant properties (k) in the theoretical equations, while the numerical model employs variable 
properties. From this comparison study, it was felt that the model was properly developed for filling and 
evacuation processes. 


Simulates Pressure Regulator by 
Pressure Gradient Input 



THERMAL 

SUBMODEL 


Figure 20: Schematic of the Gaseous System SINDA/FLUINT Model 


244 




Space Station Freedom Airlock 

During the life of the Space Static. .Freedom (SSF) t crews 

of Exits- V ehicular Activities (EV As), u, order;om^^th»Moo ^^"S«. lves for g* EVA. 

space, the EVA teams will perform a uaney evacuates the chamber by using 

Once these breathing exercises are completed, the crew enters tne airiocx, cv»,u 

a vacuum pump and journeys outside. 

If during the EVA (or some other time) , 

compression, a heat exchanger is included as part of the airlock hardware. 

As evident from equations [1] and [2], the g«es witju. mSSe 

in temperature. In addition to fluid temperature vana w J5L response \o understand how these 

JSSrS dre airkwkga^ temperature and pressum, a detailed model including all the 
aforementioned effects was developed. 

^^^R^e^l'Thow^^her^^ oftiK aMockwhichisdepic^^^r^t^ulMghnder.^22ft^ 

(2.00 meters) in diameter and 8.22 ft (2.50 ^chamber during hyperbaric operations, 
and cross-flow heat exchanger assembly wh { ch which is use d to depressurize the chamber for 

Also included in this package is a centrifugal vacuum pump wtuen is useo to ucp 

EVAs. 

bodttbe^^^ 

1) The gases within the airlock may be considered ^^^n^l^miity, indicating ideal gas 

2) At any given time, the gas within the chamber is at a uniform temperature, i nai is, me 

temperature gradients across the airlock gas. ., low velocity and laminar, convection 

3) Since .he flow inside the airlock and pmduces 

effects will be ignored. On the other hand, dunng hype pe coefficient will be the same as 

substantial flow velocities. To account for this convection, the heat transter coemcicm 

that of the Space Shuttle. ~„. t . 

A) Due to the lack of gravity, there are no natural convection enects. 

5) The gas is dry so heat transfer effects due to condensation are ignored. 

6) The gas within the airlock is not a participating media. ^ trifV ,_ f1in > 

7) The Sock is constructed out of 6061 aluminum and at a uniform temperature. 

8) The heat exchanger is modeled using the NTU method (ret. /)■ 

9) There are no heat losses to the environment 

Retire 22 shows a schematic of the SINDA/FLUINT airlock model, which is Aesanre^tite 

model show in Figure 20 with several minor the^irlock’s metal 

option and is tied to a single thermal diffusion (ome-dependenjO nod ^ts either the hyperbaric 

mass so that convective heat transfer effects can be included. The PLEN represents eunc yy* 
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charging tanks or the vacuum of space. The MFRSET has been replaced by a VFRSET (Volume Row 
Rate SET) to provide a more accurate representation of the vacuum pump. 

The hyperbaric heat exchanger is modeled with its own loop, and more specifically with a pair of 
junctions. To determine its heat removal rate, the following approach was taken. First, the NTU method 
(ref. 7,8) was applied using the characteristics of a preliminary heat exchanger design (ref. 9) and its outlet 
air temperature was determined. With the outlet thermodynamic state determined, die CHGLMP (CHanGe 
LuMP) option was used to alter the current state to the new and more accurate condition and the HTRLMP 
(HeaTeR LuMP) option was used to hold the current state. The HTRLMP option maintains the desired 
thermodynamic state by supplying an appropriate heat load at die downstream junction. 



8.22 ft 




Simulates Pressure Regulator by 



THERMAL 

SUBMODEL 


Figure 22: Schematic of the SINDA/FLUINT Airlock Model. 
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0 

pppressurization Results th^ thermal resDonse of the airlock gas as it is depressurized 

^^^i^ernodel wasfirs t used to pdf -Jf 1 LJfgX* three depressurization times were 

from 10.2 psia (70.3 kPa) to 0.5 psia ganger does not operate during this 

ssssssessses^"-" 


charge air 
reach the : 
encountered 


— nQA - c; nce th crc i S no heat transicr u» w y ™ »*v>~ — — - - 

r temperature of 70 F( r 'tic «r not in While these low temperatures may be 

-ft £SC shouldremain unaffected. 

predcited qirlock cos temperature response 


ODF^OCK IB.\ 

• OR-OCK IQ. 2 

• DR.OCK IB. 3 


9 rtlMJTC ICWIB 

?.5 HIKUTE 0CFRCS9 
m rnNun: ocpbcbs 



Figure 23: Predicted Airlock Depressurization Gas Temperature Response 


“‘“‘‘t^was next us* to dettrmim : the 

and ^ 

S"" SSSi SS SSS?- convective he« ™sfer. 

... f ^1 


Figure 24 shows the predicted remperature response <****^*“£ 
examined^These results show that heat transfer effects cm mt^nt iMyre^ * 50 °F 

%f^T42S25 ffilrtSSK S— temperature another 50 «F (28 K) 

to a comfortable range. 
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PREDICTED AIRLOCK GAS TEMPERATURE RESPONSE 
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Figure 24: Predicted Temperature Response for Three Different Heat Transfer Cases 



During normal operations, the charge air tanks will be heated and cooled by radiative heat transfer 
with the environment. It is estimated that due to this heat transfer, die temperature of the charge air may 
vary from 70 °F (294.3 K) to 0 °F (255.4 K). As indicated by equation [2], for a filling process, the final 
temperature is influenced by the incoming gas temperature, so it follows that charge air temperature will be 
important during hyperbaric operations. To examine the effect of charge air temperature on system 
response, four inlet air temperature were used: 10 °F (261 K), 30 °F (272.1 K), 50 °F (283 K), and 70 °F 
(294.3 K). The charge air was used to pressurize the chamber to 2.8 atm at 13.2 psi/min while the heat 
exchanger was operating. 

Figure 25 presents the thermal response of the airlock gas during pressurization for the four charge 
air temperatures. From these results it is clear that the temperature of the incoming air plays an important 
role in chamber temperature response. 

Airlock Model Summary and Conclusions 

From this study, it was found that convective heat transfer is important in determining airlock gas 
temperature. By using the heat exchanger , the severity of the high temperatures associated with the 
compression processes will be lessened. Finally, depress temperatures are low but should be no problem 
for suited astronauts. 


The Simplified Aid for EVA Rescue Device 

During the construction phase of the Space Station, astronauts will be required to perform many 
long-duration EVAs. While conducting these EVAs, it is possible that an astronaut may become separated 


248 



from the SSF tethers and float free. Since the Space Shuttle cannot m*™"'"** a ta^S^d* 0 " 
environment and retrieve the wayward astronaut, an alternative rescue approach must be e p y 

PREDICTED AIRLOCK COS TEMPERATURE RESPONSE 


o on . ock ib.i 
• OPLOCK ia.z 


VW' IN.CT CA3 
9flF IN-ET COS 


* BfVOCK IB 3 
■ DPLOCX IB 4 


sir xncr »s 
?BT IM-ET 005 



Figure 25: Predicted Response for Four Different Inlet Gas Temperatures 

One possible rescue method employs the Simplified Aid S 23°' 
before the EVA begins, a small thruster back pack is attache® ^^^'^^propulston 
then worn throughout the duranon of the space walk TOb tack ^kco^s^rour^ 

“ d 

guides himself to a place of safety. 

As can be shown by the governing equations for supersonic nozzle flow (ref. 1 IV the performance 
of the SAFER thrusters depends almost exclusively onthe thermodynamic Sincc j s 

working fluid at the inlet of the propulsion nozzles which is supplied from the I:*®"***. ' 

^SSSAtm the storage tank during SAFER operation, expansion effects wiU cause .die ^smre vmlun the 

tank have a significant impact on the performance, size and weight of the SAFER system. 

Devdopment ofgf Nnmenc^^od uni heat fer and 

not lend themselves to simple closed-form solutions. First, the transient cooling of the storage tanx oy 
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radiative heat transfer to deep space is described by a non-linear differential-equation, and except for a few 
idealized cases cannot be solved. The SAFER system cannot be treated as an idealized thermal radiator, 
since it also exhibits conduction and re-radiation, complicating the problem and making an analytical 
solution impossible. Second, the propulsion thrusters are supersonic nozzles and their performance is 
determined by nozzle inlet and storage tank conditions, but the storage tank pressure and temperature are 
dependent upon the mass flow rate out of the thrusters. It is clear then from the above examples, that both 
situations are quite complex and numerical models must be developed. 

Storage Tank Model 

The first situation to be examined was the radiative cooling of the storage tank gas before the 
thrusters are fired, since these results are important in the development and analysis of the thruster model. 
This model considers the heat loss of the gas and gas containment system to deep space during the EVA, 
before any propellant is used. 

Figure 26 presents a simple schematic of the propellant tank. Here, the xenon gas is contained in a 
rigid metal pressure vessel, which is protected from meteoroid impacts by a thin metal shield. The entire 
apparatus is then placed in a holding mount within the plastic (maybe fiberglass) SAFER shell and 
surrounded with Multi-Layer Insulation (MLI). The gaps between the spheres may be evacuated or filled 
with MLI. The dimensions listed on this figure are only preliminary (ref. 11) since many factors (choice 
of nozzle, amount of line heating, amount of gas), which have yet to be accurately determined, influence 
the tank size. 



SAFER 

SHELL (D = 7.5 IN) 


METEOROID 

PROTECTION (D = 7.25 IN) 
SHIELD 


PRESSURE 
VESSEL (D = 7.0 IN) 


Figure 26: Schematic of the Storage Tank 


Before the numerical model was developed, certain simplifying assumptions were made. 
Assumptions were made for both the gas and its holding vessel and these are listed below. 

1) The pressure vessel and meteoroid shield are constructed out of stainless steel. 

2) Xenon is the only gas considered. The model will be developed so that other gases can be considered. 

3) Preliminary weight estimates are; 16.0 lbm (7.25 kg) for xenon, 10.6 lbm (4.80 kg) for the pressure 
vessel, and 0.5 lbm (0.23 kg) for the meteoroid shield (ref. 11). 

4) The pressure vessel, meteoroid shield and SAFER shell are treated as concentric spheres. The 
conduction and radiation conductors for this situation are found by the method outlined in reference 9. 
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5) Only the worst case cooling situation will be con f s J^® r |^r^ t S y S s t ^ re * n ° mCldcnt 
Sun, SSF, Space Shuttle or Earth on any portion of the SAFER system. 

6) MLI is* used as the insulating material. 

?, .he quiescent fluid before thruster 

^Radiating surfaces have an effective emissivity of 0.05. 

The development of a SINDA/FLUINT caprotance is 

wv - »'™> BCK — (rcf - » is used ,0 

determine the transient cooling process. 


5mB gg ^ wha.fac.ors M-»* '£&j3E*£ 

suggest appropnate conwd methods^s^©^ cwisidered radiation as the sole means of heat transport to the 
cases were examined. Four of the five cases c the containment spheres occurs only by 

deep space environment. That is, the heat t^sfer Je considered are: xenon only (Case 

"^"effective emissivity, e‘, of 0.05. For al, oa-.hca.~ahy rejected .o the co.d 
environment of space by radiation, over a maximum period of 16 hours. 

Figure 27 presents the transient thermal ^‘^^^are dso s^ nmwizwl in Table 

S S^sS^ tial,y 

.hi pressure vessel. 

yhe Storaee Tank/Thry$K r Model f . thruster system. Here, gas from the storage tank 

Figure 28 presents a simple schematic ot the thruster sy s t u. nozz i c inlet over the 

flows iiuoa pressure regulator to ensure that a 

entire time of thruster operation. The working fluid , tn umvide thrust for the back pack, 

nozzles where it is expanded isentropicaUy m t a t ^^^ n n ^ x ^^ 0 ^g^ssures, temperatures. 
Since a thruster design (manufacturer) has yet to be chosen, the exact wonang P 

nozzle areas and flow rates of this system are not known. 

Before the Um.s.er model was developed, |he 

1) All the candidate gases (argon xenon (W roximarely 500 • 

^re^.!SrS“ S £! Arotdingdevicefendialpy is consent), U.os for ideal gases there is 
no temperature drop. 
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TEMPERATURE RESPONSE OF XENON GOS 


oncwi 1 1 CflOCi KTM EWH 4.4 CflK4 

• TVCRM 12 COSE 2 1.5 CASES 



Figure 27: Predicted Cooling of Storage Tank by Radiation for Five Insulating Cases 


Table 5: Predicted Gas Storage Temperatures at Various Times for Different Insulating 
Conditions with the Initial Temperature Set at 70 °F and the Deep Space Environment 
Set at -460 °F. 


Insulating 

Condition 

Temperature 
After 6 Hours 
(°F) 

Case 1 

0 

Case 2 

50 

Case 3 

65 

Case 4 

66 

Case 5 

68 


Temperature 
After 10 Hours 

CE1 

-25 

40 

60 

64 

67 


3) The convection in the storage tank due to mass removal may be modeled as described in reference 12. 

4) The thermal model developed to predict environmental heat leak may be incorporated into this model. 

5) The thrusters and regulator performance will be based on die most recently available data (ref 11). The 
model will be built in such a way as to accommodate changes in thruster design. 

6) The storage tank volumes of the three gases in cubic feet (meters) are; 0.223 (0.00631), 0.122 
(0.00345), 0.307 (0.00869) for argon, xenon and nitrogen, respectively. 

Figure 20 was used to model the SAFER system. The storage tank is represented by the TANK 
option which allows calculation of transient pressure and temperature changes as mass is removed from 
the tank. To account for the environmental heat loss, the TANK is tied thermally to the previously 
developed SINDA model of the storage container by the use of a convection heat transfer (HTU) tie. The 
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value of this tie is determined by the and a junction 

exhausting spherical pressure vessels. Th gu P . em nioved in the main program to alter the 

BSiSSfSSSi"* «• * *™< *» *• 


m = At*P 0 




[3] 





Figure 28: Schematic of the SAFER System. 


Ba J tg£g^ ^ «« "^ F Sy>^^"on, 

without heating of the incoming gas to examine system e P® • . ^eir individual effects on 

Sion of thTslorage tank was 8000 ps.a (55.2 MPa)and 70 F (294.3 K). 

The model was next fun to examine the flSdf^Sia^arfihout 

Figure 29 shows the predicted thruster mass flmw nte . Ae noz2 j c inlet temperature is 
heating during a one minute Mention PjnoiFsr ^'“fSmteited situation, the inlet nozzle 
maintained at the initial conditions of 70 F (294.3 K), whue it is that inlet heaung is 

would be larger and bulkier. 

To further examine what factors influence system u C ’ m70 °F 

different working fluids for identical operating conditions minme of thruster 

(294.3 K). Figure 30 presents the thermal response of the storage ^ w0 eascs Since 
operation. Here, nitrogen shows a su bs tantially [1], it follows that gases with 
the thermal response of the gas within die tank 1S ^ denressurization Since argon and xenon 

asst srasra-r »* — drop ^ 

quickly than nitrogen. 
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PREDICTED FLDUI RATES 


° venom iB. i unman »caritic 

• >CNON 1B.Z WHH HEnriNG 



Figure 29: Predicted Mass Flowrates for Heated and Unheated Inlet Gas 

PREDICTED STORAGE TANK GAS TEMPERATURE 

o\22.\ * 108.2 *100 3 



Figure 30: Predicted Thermal Response of the Storage Tank for Three Candidate Gases 
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White the results of Figure 30 sug fS J*'^ ra 7 Q?p" ( 2943 1 K)^^cun?rerapCTature), attention 
(smallest temperature difference between the iiunaljro F least amount of 

least amount of heating. 


PREDICTED NOZZLE INLET HEPTING 



Figure 


31: Predicted Heat Addition to Maintain 


the Inlet at 70°F for Three Gases 


SAFER Summary and Conclusions 

Two SINDA/FLUINT models have been n tool. Since 

hydrodynamic response of various ^mponents of *e ^ they can be easily 

the exact thruster design has yet to be chosen*e ^ek have^^m so ^ 

modified to handle any changes in working fluid, thruster aes gn, 

regulators. 

One model predicts the thermal response of ^ te*Snhmd heat teak 

K environment durrng the EVA time (sixteen hours or less). 



the required conditions. 
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CONCLUSION 


Because of its versatility, SINDA/FLUINT has been applied to a wide variety of aerospace 
problems at JSC and elsewhere. For different problems, different features of the code have been used. 
Some aspects of the code which have been utilized at JSC include transient simulation, boiling and 
condensation, two-phase flow pressure drop, slip flow, multiple submodels ^ depressunzaoon. 

Several of these applications have been described in this paper. Use of S IND A/FLUINT by the aerospace 
community and others continues to increase as the capabilities of the code expand (ref. 13). 
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